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ABSTRACT 

Context. We present a theoretical review of the effect of impact polarization of a hydrogen Htt line due to an expected proton beam 
bombardment in solar flares. 

Aims. Several observations indicate the presence of the linear polarization of the hydrogen Ha line observed near the solar limb 
above 5 % and preferentially in the radial direction. We theoretically review the problem of deceleration of the beam originating in 
the coronal reconnection site due to its interaction with the chromospheric plasma, and describe the formalism of the density matrix 
used in our description of the atomic processes and the treatment of collisional rates. 

Methods. We solve the self-consistent NLTE radiation transfer problem for the particular semiempirical chromosphere models for 
both intensity and linear polarization components of the radiation field. 

Results. In contrast to recent calculations, our results show that the energy distribution of the proton beam at Ha formation levels 
and depolarizing collisions by background electrons and protons cause a significant reduction of the effect below 0.1 %. he radiation 
transfer solution shows that tangential resonance-scattering polarization dominates over the impact polarization effect in all considered 
models. 

Conclusions. In the models studied, proton beams are unlikely to be a satisfying explanation for the observed linear polarization of 
the Ha line. 
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1. Introduction 

Observations of solar flares in the hard X-ray s pectral regio n 
indicate a non-thermal origin of this radiation dFroslI Il969h . 
There are several mechanis ms that can be i dentified as a pos¬ 
sible source of this emission (lKorchaklll967l) . Presently, there is 
wide consensus among solar physicists that the most likely ex¬ 
planation is the bremsstrahlung radiation of electron beams with 
energies of the order 10-100 keV, which are injected into denser 
layers of the solar atmosphere from a coronal reconnection site. 
The differential energy spectrum of these beams at the injection 
site is usually assumed to have a pow er-law distribution ~ 
with 6 between 3 and 5 dBrow 3ES1- 

Although there is also observational evidence of high ener¬ 
getic protons (above 10 MeV) in the y-ray spectrum, radiation 
induced by protons at energies below 1 MeV is difficult to detect 
because they do not radiate efficiently at such energies. Several 
processes that lead to electron acceleration at t he reconnection 
site c an also lead to acceleration of protons (see lOrrall & Zirk^ 
Il976l and references therein) The existence of low energy proton 
beams in solar flares is still uncertain, but it is believed that they 
may also play a significant role in flare physics. Furthermore, 
energetic electron beams cannot be us ed as a satisfactory e xpla- 
nation for all flare observations (e.g.. lDoschek et aT]ll996l) . For 
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a comparison of th e effects of electron and proton beams see 
iBrown etaH d 19901) . 

There are several techniques that can be used to de¬ 
tect low energy proton beams. The emission peak in the 
Ha core due to pro ton beam bombardment was pro posed by 
iHenoux et al.l dl993l) . It has been recently shown by IXu et’aP 
d2005al) that this effect does not exist. Another technique is 
based on measurements of the red-shifted emission in line 
wings, especially in the hydrogen Lya line, which is a con- 


seauence of the charse exchanse effect (lOrrall & Zirked 

19761 

Canfield & Chang 

119851; iBrosius et all 119951; iFang et all 

19951 

Zhao et al. 19981 

Brosius & Woodgatel 19991). A different an- 


proach is based on the fact that anisotropic excitation of the 
chromosphere atoms by a directed proton (or electron) beam in¬ 
duces a preferential population of particular Zeeman sublevels, 
i.e., the impact atomic polarization and, consequently, an emis¬ 
sion of linearly polarized radiation dHenoux et al.ll 199^ . 


Some observations indicate the existence of linear polariza¬ 
tion of the Ha line in solar flares above 5 % or even as high as 


10 %( 

Henoux & Chambell 99011 Vogt & Henouxll996l;IXu et al.l 

12005b 

) and preferentially oriented towards the disk center (it 


will be denoted as radial in this paper) and also parallel to 
the limb (i.e., tangential). This effect is usually interpreted as 
a consequence of the impact polarization by a vertical proton 
beam with energy below 200 keV in the Ha core forming layers. 
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Contrary to these measurements, there are observati ons that indi¬ 
cate n o linear polarization in a wide range of flares (iBianda et alJ 
l2005h and argue for isotropization of the beam at crucial at¬ 
mospheric depths. Studies of Ha polarization in solar promi¬ 
nences show that transitions between fine structure levels within 
the shell n = 3 of hydrogen caused by collisions with back¬ 
ground electtons_andjprotons lead to a significant depolarizing 
effect (iBommier et al.lll986bh even at densities of the order of 
10'° cm“^. The electron (proton) density in a flaring solar chro¬ 
mosphere at the Ha formation levels is above lO'^cm^^, there¬ 
fore this effect could significantly reduce measurable polariza¬ 
tion. 

In recent years, a first quantitative estimation of the hydrogen 
impact polarizatio n un der the flar e cond itions has been done by 
IVogt et aP (Il997h and IVoet et alJ (1200 ll) . In these calculations, 
a self-consistent radiative transfer code for unpolarized radia¬ 
tion has been used to find the hydrogen ionization degree and 
the radiation intensity at unit optical depth of Ha for several 
semiempirical atmospheric models and under the effect of pro¬ 
ton beams having several parameters. In these calculations, the 
following approximation has been used: The equations of statis¬ 
tical equilibrium and the radiation transfer equations have been 
decoupled and the Ha line has been assumed to be optically thin. 
The proton beam energy distribution at an injection site has been 
assumed to be similar to the one deduced for electron beams, 
i.e., a power law. The lower energy cut-off Ec is usually set be¬ 
tween 100 to 200 keV. This energy approximately corresponds 
to the energy necessary for a proton to reach the Ha line-forming 
layers. Using an assumption that the power-law distribution pre¬ 
serves its character while the beam propagates to the upper chro¬ 
mosphere, these authors calculated the polarization degree of 
emitted radiation. Depolarization by collisions with electrons 
and protons of the ambient medium have also been taken into 
account. Depending on the atmosphere model and beam param¬ 
eters, the predicted polarization degree has been up to 4.5 %. 
However, a complete solution of the coupling of the polarized 
radiative transfer to the atomic equilibrium has not been tackled. 

The aim of this paper is to present a theoretical analysis of 
the proton impact polarization phenomenon in the hydrogen Ha 
line based on NLTE polarized radiative transfer. Our goal is to 
verify the assumption that the measured linear polarization can 
be interpreted as due to this mechanism. In particular, we have 
focused our attention to a limiting case of unidirectional non¬ 
deflecting beams, which are supposed to generate highest polar¬ 
ization due to their extremal anisotropy. For simplicity we use a 
static ID model of the flaring chromosphere. 

In the second section, the deceleration effect of the chromo¬ 
sphere on proton beams is reviewed and some conclusions are 
made about the low-energy beam fluxes in Ha-forming layers. 
These are compared with the assumptions made in the previ¬ 
ous works. The third section describes our NLTE solution of the 
unpolarized transfer in an atmosphere affected by non-thermal 
excitation an d ionization by proton beam. This solution is simi¬ 
lar to that of iKasparova & Heinzell (l2002ll . which has been per¬ 
formed for electron beams. We use the same approach to cal¬ 
culate the volume densities of thermal electrons and protons in 
the chromosphere. These densities are then used in our polar¬ 
ized transfer code. After that we describe the framework of the 
quantum density matrix and equations of atomic statistical equi¬ 
librium on the basis of irreducible tensorial operators. The fifth 
section is dedicated to the problem of collisional rates for all 
transitions used in our modeling. S ect, [^ contains a comparison 
of results obtained by IVogt et aH (1200 ih with our calculations 
in the last scattering approximation. A brief description of po¬ 


larized transfer solution and the method used in our calculations 
can be found in Sect.|7] The results and their discussion are sum¬ 
marized in Sect. 0 In the last section our conclusions are made 
and the main consequences for the interpretation of observations 
are pointed out. All expressions in this paper are written using 
the CGS system of units. 

2. Solar atmosphere models and beam propagation 

2.1. Proton beam propagation 

Let us briefly review the interaction of a proton beam with the 
solar chromosphere. We assume a beam in the coronal reconnec¬ 
tion site, which vertically penetrates the chromosphere. A hori¬ 
zontal motion of the protons is neglected to obtain as anisotropic 
a velocity distribution as possible. The energy of protons neces¬ 
sary to reach the Ha formation levels is of the order of 100 keV 
at the top of the chromosphere (see below). Protons of such 
high energy do not lose a lot o f energy in the inter action with 
extremely hot coronal plasma (iHenoux et al.l 1 199311 . hence we 
may neglect any interaction with the coronal mass in this paper. 
In all our calculations, the energy of the superthermal protons 
(> 1 keV) of the beam is high above the average energy of ther¬ 
mal motions in the chromosphere where Ha is formed, thus w e 
can use the so-called cold target approximation (lEmsli ell 19781) . 

In a partially ionized medium, a superthermal proton beam is 
decelerated by collisions with the background electrons and ions 
and also by elastic and inelastic collisions with neutral atoms, 
especially the atmos phere’s r nain c onstituent, hydrogen. We will 
use the approach of lEmslid (Il978l) to describe this deceleration 
and the beam energy deposition into the atmosphere. The decel¬ 
eration by charged electrons and protons can be quantified by 
means of the Coulomb logarithm A. The inelastic and elastic 
scattering on the neutral hydrogen is similarly described by A' 
and A", respectively. These logarithms vary slightly according 
to atmospheric properties, but these changes are small within a 
wide range of physical characteristics and it will not affect our 
results significantly if we suppose these quantities to be con¬ 
stant. The values adopted in this paper have been set according to 
typical physical properties of the upper chromosphere: A = 23, 
A' = 3. Elastic scattering of the proton beam on neutral hydro¬ 
gen is negligible in comparison to other processes and will be 
neglected. 

Let F(E,N) be the energy distribution of the number flux 
of beam particles at the column depth N. At the injection 
level we set N - 0. Th e general form of this distribution is 
(ICanfield & ChandfT^ 

F{E,N) ^ F{Je^ + El,0)—^= . (1) 

For the purposes of our calculations, we have used an ini¬ 
tial power-law energy distribution at the injection site, which 
is usually con sidered in the non-thermal flare heating problems 
(iBrownll 19711) 

F{E, 0) = (5 - 2)El-^£,QE-^e{E - E^) , (2) 

where 

is the Heaviside function, Ec is the lower energy cut-off of the 
distribution, and So is the total initial energy flux of the beam. 
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Ec is usually assumed to lie between 100 keV and 200 keV, ac¬ 
cording to minimal energy necessary for protons to reach the Ha 
formation layers. The formula dB for th e particular case of initial 
distribution (|2]) gives (IZhao et Sll998h 


F(E, N)^ (6- 2)Et^6o ^ ^ , for £ > , (4) 

(£2 + £2 )- 

where 

E^(N) = I < £c) (5) 

[ 0, {En > £c) 

is the minimum energy of the distribution at a given depth. 

The energy deposition rate into the neutral hydrogen is given 
by 


Im{N) 


Kny 


■(b-2)(l -.x)A' 


tWp £() 

tMe Er 




where 




N/N„ (N<N,) 
1, (N > Ac) 


( 6 ) 

(7) 


and Bx(a, b) is the incomplete beta function. Ac is the depth that 
can be reached by protons w ith the initial energy £c. Eg. (Ibl ) is a 
special case of the Eq. (1) of lKasnarova & Heinzell (l2002h . here 
for a non-deflecting vertical proton beam. This energy deposition 
rate has been used in calculations of the ionization degree of the 
chromosphere described in Sect. [3] 


2.2. Flux distribution at Ha formation level 

The distribution (IH has been plotted in Fig.[T]for several values 
of En, typical energy cut-off £c = 150 keV and two initial spec¬ 
tral indices 5. The first fact that should be pointed out is that the 
energy cut-off does not significantly depend on E^ if E^ < £c 
until En ss £c. It is because high energy protons are decelerated 
inefficiently. Only if £^ ^ £c, protons are slowed to energies 
close to 5 keVQ- Using Eq. (|4]i we And that the flux of the beam 
at the local energy cut-off E^ is lower by factor {EcIEmYE^IEn 
compared to the flux £(£c, 0). Hence we have a strong decrease 
of flux close to the interesting energy range (i.e., E^ <si E^) by 
the factor of approximately E^jE^. The total flux is rapidly re¬ 
duced in deeper layers > £c). It can be shown that the flux 
maximum is at the energy Ef^j V^. After crossing some critical 
depth, the flux of the beam with lower 6 dominates the one with 
higher 6 at all energies - because initially there is a higher num¬ 
ber of high energy protons in the small-b case, which are now 
decelerated to low energies. As a result, energy is more effec¬ 
tively deposited by small-i5 beams in lower depths, while high-5 
beams are decelerated in upper layers. In any case, the decrease 
of flux dF(E,Eff)fdEN is steep for any 5 and one could expect 
that impact polarization will be sensitive to the E^ value because 
there is only a small A£jv interval in which low-energy flux is 
not negligible, and this depth interval should overlap with the 
Ha line center formation region as much as possible. It is not 
possible to have a power-law-like distribution at this layer with 
a local cut-off of about 5 keV for beams starting at the top of 
chromosphere with energies above 100 keV or higher. 

This leads us to the second conclusion. The energy distri¬ 
bution of the flux cannot be approximated by the power-law 

* As we will show later, the low-energy protons (4-5 keV) are most 
effective in producing impact polarization of « = 3 hydrogen level. 



Fig. 1. Energy distribution of the beam with initial energy cut¬ 
off £c = 150 keV at depths with different Solid curves corre¬ 
spond to spectral index 5 = 5, dashed curves to 5 = 3. Erom up¬ 
per to lower, the curves are plotted for £a^ = 0, 145 and 230 keV. 
The cut-off energy at depths close to Ac is still of the order of 
£c. At deeper depths the distribution 5 = 3 dominates over 5 = 5 
at all energies. 


curve if Eat is close to Ec. If these energies are comparable (and 
that is the case of our interest here because we have for max¬ 
imum impact polarization E^ ~ 5 keV <si Ec), the use of the 
power-law distribution leads to an unrealistic overestimation of 
th e polarization deg ree. Su ch an approxima tion has been used 
bv IVogt et al.l (1 19971) and bv IVogt et al.l (1200 Ih with lower energy 
cut-off varied from 1 keV to 20 keV and with a total beam num¬ 
ber flux and 5 conserved. The theoretical values of polarization 
obtained from these calculations should therefore be revised. 


3. Unpolarized radiation modeling and background 
electron and proton densities 

Our solution of unpolariz e d rad iation transfer follows the 
approa ch of iHenoux et al.l (1 19931) and iKasnarova & Heinzell 
(l2002l) . We take a ID static atmosphere with fixed tem¬ 
perature structure and other plasma properties given by the 
sem iempirical chromosp here model of t he flaring atmosp here 
El (iMachado et alll 19801) . and following IVogt et S ] (11997 1) and 
Vogt et al. ( 20011) we also use the model VAL E ( Vernazza et al.l 


19811) for comparison. These models have been created using a 


number of line and continuum observations. 

The non-thermal proton beam dissipates its energy while 
propagating through the matter of the chromosphere causing 
its heating and the modification of atomic level populations 
and ionization. The non-thermal heating is already included in 
the semiempirical flare models to explain the observed emis- 
sion. Our approach to t he modeling follows t he ap proach of 
iHenoux et al.l (1 19931) and lKasnarova & Heinzell (1200^ . We con¬ 
sider a fixed temperature structure of the atmosphere as given 
by the semiempirical atmosphere models. The pressure and sta¬ 
tistical equilibrium equations are solved together with radiation 
transfer. Then we study an influence of the non-thermal colli- 
sional rates on the line profiles in comparison to the thermal 
model (differential approach). It is important to notice that the 
temperature of the flaring atmosphere described by semiempir¬ 
ical models is in general overestimated because it has been de¬ 
termined to explain the increased radiation emission regardless 
of the non-thermal excitation and ionization. The same approach 
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has been used for the solution of the polarized radiative transfer 
(Sect.|7]i to obtain the changes of the Ha ^//-profiles after intro¬ 
duction of the non-thermal collisional rates into statistical equi¬ 
librium equations. The only difference is that the polarized solu¬ 
tion uses the electron and proton densities precalculated by the 
unpolarized transfer code. W e have used the same MALI code as 
iKasparova & Heinzell (l2002h . now adapted to treat proton beam 
bombardment described in Sect. |3 The code has been run with 
a four-level plus continuum hydrogen atomic model and non- 
thermal collisional rates have been calculated from the energy 
deposi tion rate into hydroge n (|6]l using the expressions (10) and 
(11) of lHenoux et al.l (|l993l) . 

The electron densities (which are close to proton densities at 
Ha core formation layers) are plotted in Figs. |2] and [3] for both 
models under consideration, several beam fluxes, and three dif¬ 
ferent values of 6. These calcul ations show that am bient electron 
(and proton) densities used by IVogt et aP (Il997h are underesti¬ 
mated (cf. Tables 2 and 3 therein) and the same is true for the 
mean radiation intensities. These values have been computed in 
the layers of total (coronal + chromospheric above Ha-forming 
layers ) column mass depth 2 .484x 10~^ g cm“^ (VAL F, see Table 
15 of I Vernaz za et al.l 1 198 Ih an d 3.186 x 10 gem ^ (FI, see 
Table 3 of Machado et alJl9^ . where the optical depth of Ha is 
approximately unity. Because our backgroun d particle densitie s 
did not correspond to the ones presented by IVoet et akl (Il997h . 
we have used their radiation transfer code which is based on the 
standard A-iteration process. We have found that this code set¬ 
ting (the accuracy of 10“^ in maximum relative change of atomic 
populations between subsequent iterations) leads to insufficient 
convergence and consequently to the underestimation of elec¬ 
tron and proton densities. An increase of accuracy led to values 
close to ours, but at the cost of an extremely large number of A - 
iterations (for similar tests see also lKasnarova & Heinzel|[2002h . 

In Fig. |4] one can see the theoretical Ha line profile for the 
same beam fl uxes as in Fig.lH T hese profiles are consistent with 
the result s of IXu et al.l (l2005ah . In comparison to profiles pre¬ 
sented bv iHenoux et al.Nl993l) . the central emission in the line 
core does not appear. 


4. Formalism of density matrix and local 
equilibrium 

The electron densities in the layers of interest do n ot ex- 
ceed 10'^ cm“^. F ollowing iBommier & Sahal-Brecholj (Il982h 
and [Bommier et all (Il986al) the hyperfine splitting of the hydro¬ 
gen levels is completely negligible. In fact the hyperfine splitting 
of the relatively long-living 3s^^^ level (about 50 MHz) is higher 
than its inverse radiative lifetime (of the order of 1 MHz) and 
thus th e hyperfine leve l s do no t overlap. However, it was veri¬ 
fied by iBommier et alJ (Il986ah that taking the hyperfine struc¬ 
ture into account does not affect the linear polarization of the 
Balmer lines in any significant way. Moreover, in the conditions 
under consideration, the lifetimes of the levels are strongly re¬ 
duced by collisions with the charged perturbers. The lifetimes of 
the n - 3 levels are of the order of a few 100 MHz in the Ha 
cor e-forming layers due to the dipolar transitions nlj nl ± 1 f 
(cf. lSahal-Brechot et alJll9^ : hence, much greater than hyper¬ 
fine splitting. Therefore, we take into account only the fine struc¬ 
ture splitting of the levels. We suppose that there are no quantum 
coherences between different energy levels. This assumption is 
valid in our model of the hydrogen atom where the individual 
level widths are smaller than distances between them or where 
the selection rules for dipolar optical transitions prevent creation 



0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 
M [g cm“^] 



0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 
M [g cm“^] 



0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 
M [g cm“^] 

Fig. 2. Electron density for the model FI. The plots correspond 
to the case 6-3 (upper panel), 6-4 (middle panel), and <5 = 5 
(lower panel). The thermal case is plotted by a solid line, the 
non-thermal beam fluxes are So — 10^ (dot), 10^ (dash), 10'° 
(dash-dot), and lO" ergem^^ s“' (dash-dots). All beams have 
= 150keV. 

of coherences. As a result, the fine structure levels are supposed 
to be completely separated. 
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Fig. 3. Same as Fig.|2j here for the VAL F model. 


To describe the state of hydrogen atoms in an equilibrium 
with the radiation field and charged colliders , we adopt the 
framework of the atomic density matrix p (e.g., lFanolll957h . In 
the dyadic basis \nl jm) of hydrogen states, we adopt the common 
notation: n is the principal quantum number, / the orbital one, 
and j is the total angular momentum. Magnetic quantum number 


5x10 



IxlO b_^^^^^^^^^^^^^^^_d 

-4-20 2 4 

AX [A] 

Fig. 4. Ha line profiles for the model FI and S - A. The partic¬ 
ular profiles correspond to the same beam fluxes as in Fig.|2l 


is represented by m. The natural basis for the density matrix op¬ 
erator in polarization studies is that of irre ducible tensorial oper - 
ators (ITO or the multipole expansion) (ISahal-Brechodl97^ . 
In all our development we suppose that all the particle velocity 
distributions and the radiation field are axially symmetric with 
the symmetry axis in a vertical direction. We also assume that 
the magnetic field lines are oriented along the vertical axis and 
that the strength of the magnetic field is at most a few hundred 
gauss (thus we can neglect the Zeeman splitting of the levels). 
All of these assumptions lead to a great simplification of calcula¬ 
tions because most of the density matrix elements vanish: All the 
coherences between the wave functions of the states vanish. In 
the basis of ITO only the components remain. These diago¬ 
nal elements represent population (k = 0) and alignment (k - 2) 
components of the different levels. Again, due to symmetry rea¬ 
sons (pnijm - Pnij-m), the elements with odd k v anish. One could 
take into account the anti-level-crossing effect (lBommieiill980h 
related to the so-called align ment-to-orientation mechanism (cf. 
iLandi DeglTnnocentill982h . which can lead to creation of some 
non-zero density matrix components of odd rank (so-called ori¬ 
entation), even if no circularly polarized radiation is present. In 
fact, this effect is of little importance in our current model be¬ 
cause the magnitude of the created orientation remains well be¬ 
low the magnitude of the alignment. Therefore we can safely 
exclude this effect from our study. The higher multipoles {k > 2) 
do not affect the solution under given the physical conditions 
and they have been neglected^ We limit our analysis to the three 
principals levels of the hydrogen atom. 

Our calculations concern a static model. The local hydrogen 
equilibrium is expressed by the equations of statistical equilib¬ 
rium (ESE) 

Hp = M . (8) 

The elements of matrix 11 consist (in the impact approximation) 
of the collisional and radiative rate s simply summed together: 

. = . H- . dBommier & Sahal-Brechotl 

Il99lh . The diagonal elements in nlj = n'l'f stand for relax¬ 
ation of the level nl j to another levels, while the nondiagonal el¬ 
ements stand for populating a given multipole component. The 


^ The only non-zero multipole A; > 4 in the third-principal level hy¬ 
drogen is ^‘'^Pq and one has 
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relaxation rates will be marked by the index “R” in this pa- 
pe r. The detailed analy sis of the structure of 11 can be found 
in LSahal-Brechol (Il977h . p is the formal vector of density ma¬ 
trix components "’■'Pq. The normalization condition on the den¬ 
sity matrix has to be introduced; The sum of all level populations 
is equal to 1. We know that the relative population of level nlj is 
equal to and hence we may replace the first row of 

n by the appropriate constant elements and the right-hand-side 
of Eq. (I8]l reads u = (1,0,..., 0)^. 

The particular expressions for the R-matrix elements can be 
found in lLandi DeglTnnocentil(ll984ll . They are calculated in the 
lowest order of quantum electrodynamics from a known radia¬ 
tion field. Only the single events of emission and absorption are 
considered. We adopt the approximation that there is no coher¬ 
ence between the frequency of absorbed and emitted photons. 
This approximation is known as a complete frequency redistri¬ 
bution (CRD; cf. also Sect. |2ll. The elements of the matrix of 
collisional rates C can be calculated taking into account several 
collisional processes as will be described in the following sec¬ 
tion. 


5. Collisional rates in detail 

First of all we stress that we limit our analysis of collisional tran¬ 
sitions to density matrix multipoles of k = 0 and k - 2. Higher 
ranks do not significantly affect polarization in the studied 
problem and will be neglected. Collisional rates and 

^nJjtn'i'f hydrogen can be calculated from cross-sections 
of appropriate transitions by integration over the 
relative velocity distribution /(v) of particles 

-JVpJ , (9) 

and similarly for relaxation rates. A^p stands for the perturbers 
volume density. 

The processes taken into account in our analysis are the fol¬ 
lowing; 

1. Fine structure dipole transitions induced by ambient elec¬ 
trons g- -I- (Hi)”'^' e- + 

2. Fine structure dipole transitions induced by ambient protons 

-H (Hi)"'-'' ^ -h (Hi)”'*'^'. 

3. (De)excitation by ambient electrons + (Hi)"^-' + 

(H I)"'''-''. 

4. Excitation by proton beam pg -i- H i —» pg -H (H i)*. 

5. Charge exchange excitation pg + H i ^ (H i)* H- p^. 

The index “B” is reserved for particles of the beam and an as¬ 
terisk (*) is used for an excited state of the neutral hydrogen 
H I. p^ and denote protons and electrons, respectively. The 
other possible collisional transitions have been neglected due to 
their negligible effect. The cross-sections of these processes have 
been obtained in several ways. 

The fine structure transitions within the same shell {nlj —> 
nl ± I f) hav e been calculated using the semiclassical perturba¬ 
tion method (ISahal-Brechot et al.lll99^ . This approach is accu¬ 
rate for calculation of transitions between close levels (in com¬ 
parison to collision energy) in a dipolar approximation. Using 
the formalism of the previously referenced paper, rates can be 
expressed in the form 


, ( 2 ) ( 2 ) 

-I- c, ■ -.a ,■ , 


R^k^k' _ AT t ki 

^nlj^n'I’f - ^Pydkk’ 

( 2 ) 


+ r 

k—^k'j-^f nij—^n 


'I’j’] ’ 

:o) 

i/ j^n'l 


with 


O', 


(0) 

nl j—>n‘ 


X co 




.( 2 ) 

nlj-tn'l’j’ 



{An 


42) 


■7;^„,,,Xv)/2(v)vMv 


( 10 ) 


( 11 ) 


( 12 ) 

(13) 


The cross sections have been defined by expressions 

(59) and (60) of the referenced paper, /o and /2 are the monopole 
and quadrupole components of the relative velocity distribution 
with axial symmetry. The angular coefficients c have been de - 
fined by expressions (66)-(69) of ISahaTBrechot et al.l (Il996ll . 
Although this method is limited only to treatment of dipolar tran¬ 
sitions I ^ I ±1, these transitions are dominant in the important 
range of collisional energies at chromospheric temperatures. 

Excitation transition probabilities from the ground state lii 
to upper levels induced by thermal electrons have been calcu¬ 
lated using the cross-section data of database AMDIS0 These 
data are provided for Is nl transitions neglecting the fine 
structure of the levels, and therefore we use an approximation by 
the angular coeff icient to obtain the fine structure cross-sections 
(IVogt et alJl2001h . 




/ / k 

j j 2 


cr 


Is^nl ' 


(14) 


At the present time, we do not have any adequate cross-section 
data for the fine structure transitions between levels n - 2 and 
« = 3. Therefore we have calculated these cross sections us¬ 
ing the same semiclassical perturbation method, although these 
cross sections are overestimated and cannot be used for non- 
dipolar transitions 2s 3d. 

The cross sections of direct excitation to levels n - 2 
and » - 3 by protons hav e been calculated using the data of 
[Balanca & Feautri^ (Il998ll . These authors give cross sections 
(T^^ni population and alignment excitation in the energy 
range of 1 keV to 100 keV. As well as in the case of excitation by 
thermal electrons, we use the expression (fT4l) to obtain data for 
the fine structure transitions. The transitions n - 2 n - 3 
have been neglected. The self-consistent solution of radiation 
transfer requires us to know the cross sections at energies above 
100 keV. This is due to the behavior of the energy distribution 
of proton beams at different depths. Although the total flux of 
the beam decreases quickly as the depth exceeds Nc, the role of 
high energy protons increases; the maximum of the distribution 
F{E, N > Nc) is at the energy E^/ s/d. This leads to the emission 
in the near Ha line wings and neglecting these energies leads 
to unrealistic line profiles. We have used the semiclassical per¬ 
turbation method to calculate dipolar cross sections at energies 
above 100 keV. The transitions li —> 3t/ cannot be calculated 
b y this method, and therefore we used a similar approximation 
as IVoet et aP (1 19971) and set these transitions to one tenth of the 
value of Is ^ 3p cross-sections. This approximation is based 
on the observation that the electron cross-sections for 1 s ^ 3t/ 
excitation are approximately 1/10 of Is —> 3p at high energies 


^k—*k' 
^nlj^n'V f 



(0) ^(0) 
k.i^y nlj^n'l'j' 


^ Atomic and Molecular Data Information System, 
http://www-amdis.iaea.org/ 
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( Aboudarham et afll992h. and it is also true for 100 keV protons 
(iBalanca & Feautri^ll998h . While this is not fully justifiable, it 
should improve accuracy of the solution. 

The effect of the char ge exchange has also been t aken into 
account using the data of IBalanca & Feautrierl (Il998h . On the 
one hand this process is efficient at energies around 25 keV, on 
the other hand its cross section decreases fast at higher ener¬ 
gies. A more detailed analysis of the effect of the charge ex¬ 
change and the associated Doppler shift of the line can be found 
in Sect. [T] Contrary to direct excitation, no approximation of the 
cross-sections at energies above 100 keV has been used. These 
cross-sections can be safely neglected due to their small magni¬ 
tude. 


6. Impact polarization without polarized transfer 

In this section, we will c ompare several re sults of atomic impact 
polarization obtained by I Vogt et al.l (1200 ih with our calculations 
performed by a similar method. We will consider the emission of 
radiation under the proton beam bombardment in chromospheric 
conditions. The local atmosphere conditions (temperature, ther¬ 
mal electron and proton densities, and mean radiation intensity 
in the hydrogen La, LyS, and Ha lines) are obtained from the un¬ 
polarized radiative transfer solution. These co nditions are taken 
at uni t Ha line center optical depth. Following IVogt et^ (1 19971 
1200 ill we assume that proton beam energy distribution is given 
by a power law with the same spectral index S as at the injec¬ 
tion coronal site. The emitted radiation, which is polarized due 
to beam impacts, is assumed to be unaffected by radiation trans¬ 
fer in upper layers in this approach. Although a rough approxi¬ 
mation, it can serve as a comparison of our calculation of colli- 
sional rates and local equilibrium, which will be used in further 
self-consistent solutions. 

After calculation of the density matrix elements from Eq. (O 
one can calculate the ratio ep/e/ of emission coefficients of the 
Stokes parameters Q and I. When observed at solar limb (i.e., 
90 degrees from the vertical direction) and through the optically 
thin layer, this ratio could provide an estimation of the emergent 
light polarization degree. From known atomic density matrix el¬ 
ements, the polarization of emergent radiation reads 


P9o(n 

with 


«') = - 




,( 2 ) 


2 V277^“* , - , 

' n^n' n-^n' 


(15) 


§ 


(k) 

H—>H' 


^(_l)iVV'(2j + I { { I "'"Po ■ (16) 


The quantity A„ij^n'i ’/ is the Einstein coefficient for spontaneous 
emission and the entities in the braces are the so-called 6 j sym¬ 
bols. At the level of unit Ha optical depth, the beam energy cut¬ 
off is, by definition, close to ze ro. This cutoff has been varied in 
the interval 1 keV to 20 keV bv IVogt et alJ (1200 ill to find the de¬ 
pendence of the polarization degree. These authors show that the 
polarization degree can be expected as high as 2.5 % for model 
VAL F and 0.7 % for model El (6 = 4), or 5 % for model VAL F 
and 1.2% for model El (6 - 5). The results of IVogt et al.l (1200 111 
contained in their Table 1 are plotted in Fig.|5] 

We have repeated these calculations for the same conditions, 
with Ecirua = 1) varied between 1 keV and 20 keV and to¬ 
tal initial beam energy flux between 10* to 10'^ erg s“'cm“^. 
The physical conditions of the atmosphere have been taken from 



Fig. 5. Emergent polariz ation degree (in pe rcent) of the on-limb 
observation calculated bv IVogt et al.l (1200 Ih in the last scattering 
approximation. The horizontal axis shows the value of the local 
energy cut-off, and the vertical axis shows the values of the initial 
beam flux. 
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Fig. 6. Same as Fig. |5j here calculated using another numeri¬ 
cal integration technique and different cross-section data for the 
electron excitation (see text for details). 


Tables 2 and 3 of IVogt et al.l (Il997h . The corresponding polar¬ 
ization degree for VAL F model is plotted in Fig. |6] 

Although the order of both results agrees, there are differ¬ 
ences between these solutions. The most important difference 
is a steeper decrease of polarization degree at higher energies. 
The reasons for these differences are the different cross-sections 
used for electron-hydrogen excitation and especially the differ¬ 
ent technique of numerical integration used for calculation of 
non-thermal collisional rates. While the Gauss-Laguerre quadra¬ 
ture has been used in Vogt’s calculations, we have chosen an 
adaptive Simpson rule with a very fine energy grid refinement 
at low energies, wh ere fast changes of cross-s ection data play a 
role (cf. Fig. 11c of IBalanca & Feautrieiill998ll 

A further step of our analysis has been the evaluation of the 
polarization degree in the same approximation but under local 
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Fig. 7. Same as Fig. |6l here calculated with correct densities 
of the background electrons (and protons) and modified mean 
intensities of the radiation field. 


physical conditions calculated using the correct solution of the 
unpolarized radiation transfer (see Sect. [3]l. An increase of the 
mean radiation intensity and higher depolarization due to in¬ 
creased ionization degree lead to a significant decrease of the 
polarization degree by more than one order of magnitude (see 
Fig. El). The expected polarization degree is very sensitive to 
nlj nl ± 1 f depolarizing transitions between fine structure 
levels of the same shell caused by background perturbers. The 
results of this section show that the line polarization is in fact 
smaller than the previously reported calculations, even if we use 
the not fully correct approximation (EJi for the energy distribu¬ 
tion of the beam at the Ha core-forming layer. 


7. Polarized radiation transfer solution 

A self-consistent NLTE solution of the polarized transfer may 
provide additional information about the polarization degree of 
the Ha line and the effect of the proton beam. In contrast to the 
results of the previous section, it also provides complete emer¬ 
gent intensity and fractional Qjl polarization profiles. In this 
section, a short description of the polarized transfer solution is 
contained. It is beyond the scope of this paper to review the com¬ 
plex problem of self-consistent modeling of the polarized trans¬ 
fer in hydrogen lines and thus we limit our explanation to a few 
notes about the method. The detailed description of this problem 
is a subject of other papers (ISter)anll2006l 2007, in preparation). 

It is advantageous to describe the polarized radiation 
field by means of the so-called Stokes vector S(v, £2) = 

(/, Q, U, V)^(v, fl), in which / is the specific intensity, Q and U 
correspond to the linear polarization parameters with respect to 
two axes, and V stands for the circular polarization. We limit 
our analysis to the cylindrically symmetric problem where the 
natural choice of the reference frame leads to a reduction of the 
Stokes vector to only two nonvanishing parameters, S(v, il) - 
(/, Q)^(v, fl), while [/ and V are identically zero. Stokes param¬ 
eters are computed from the radiation transfer equation (RTF) 
formally identical to the unpolarized transfer case (all depen¬ 


dences on the radiation frequency v and direction of propagation 
n are suppressed): 

dS 

— =6-(i7-i/,)S. (17) 

ds 

The 2-component vector e stands for the spontaneous emis¬ 
sion in Stokes parameters, r/ is the 2x2 symmetry matrix 
of the absorption, and rjs matrix stands for the stimulated 
emission. Explicit fo r ms of these quantities can be found in 
iLandi DeglTnnocentil (Il984l) . Eor the purposes of this paper, 
it is only important to notice that values of these quantities 
can be calculated directly from the atomic density matrix co¬ 
efficients p^ij. In the solar chromosphere, nonlocal coupling of 
atomic states between different points due to radiation is strong. 
Therefore we have to solve the nonlinear system of RTE (fTTl i 
together with ESE ([8]) at every point of the atmosphere. This 
so-ca lled NLTE problem of the 2nd kind (iLandi DeglTnnocentil 
Il987h can be solved numerically using efficient iterative NLTE 
solvers. 

We have used our new multigrid technique that calculates 
the Stokes profiles for a given temperature and density struc¬ 
ture of the atmosphere for an arbitrary line overlapping, and the 
CRD approximation. This method uses the accelerated lambda 
iteration technique applied to polarized radiation together with 
the m ultigrid acceleration technique (cf. iFabiani Bendicho et al.l 
1 19971) . which improves the convergence rate from to 

0(N), where N is the number of spatial discretization points in 
the atmosphere). The formal solve r of our method is bas ed on 
the short characteristics approach (lOlson & Kunaszl[l987l) with 
the parabolic interpolation of the source function. 

The effect of stimulated emission is taken into account as 
well as the continuum opacity and emissivity. In the present 
calculations, we have restricted our solution to the first three- 
principal levels of the hydrogen atom and the transfer has been 
solved in three spectral lines, i.e.. La, L/3, and Ha. A 50-point 
quadrature in frequency for each line and a 14-point quadrature 
for ray directions has been applied. 

Finally, one simplification of our approach should be men¬ 
tioned. While in common situations the density matrix is treated 
as independent of atomic velocity, this is not the case for 
anisotropic charge exchange interaction. In this case, there is a 
systematic Dop pler shift in emission dependent on the dir ection 
of observations (ICanfield & Changll985l:lFang et al.lll995l) . This 
effect is dominant in the La and L/3 lines, while the Ha line pro¬ 
file is not much affected. In our approximation we do not take 
into account any Doppler shifts and assume charge exchange to 
be symmetric with respect to the center of the line. This approx¬ 
imation could lead to a slight overestimation of the impact po¬ 
larization effect on the Ha line. 

8. Results & discussion 

Particular results of the Ha polarization in the approximation 
neglecting polarized transfer have been presented in Sect. |6] 
Using our multigrid method, in Sect. E] we have solved the self- 
consistent NLTE transfer problem of the 2nd kind for a grid of 
atmospheres based on El and VAL F models affected by proton 
beams of different fluxes and energy distributions. 

Figs. Eland |9] show the theoretical emergent linear polariza¬ 
tion profiles Qjl of the Ha line for the FI and VAL F models 
respectively. The intensity /-profiles (not shown) are similar to 
the profiles plotted in Fig. |4l although the intensity in the line 
core and near wings is overestimated by factor of approximately 
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Fig. 8. Emergent fractional linear polarization Qjl profiles com¬ 
puted close to limb ip. — 0.11) for the FI model and differ¬ 
ent beam energy distributions: 6 = 3 (upper), 6 = A (middle), 
and 6 = 5 (lower). The different lines correspond to the same 
beam fluxes as in Fig.|2 All the calculations have been done for 
£c = 150keV. 


1.3-1.5. These differences between intensity profiles are due to 


Fig. 9. Same as Fig.[8j here for the VAL F model. 


different collisional rates used in both calculations and because 
the polarized solution was restricted to three-principal levels of 
the hydrogen atom. However, neglecting the upper level should 
not affect the emergent polarization profile much. In our nota¬ 
tion, a positive sign in 2/7 plots means the tangential direction 
of polarization, while a negative sign is the radial one. Because 
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of the symmetry reasons the highest fractional polarization is 
expected at the limb (p = 0) and the degree of polarization de¬ 
creases to zero at the disk center (// = 1). 

The resonance polarization of the Ho- line in thermal atmo¬ 
sphere (no beams) shows a small tangential polarization of de¬ 
gree 0.01 % for FI and 0.4 % for VAL F. This shows the domi¬ 
nant role of resonance polarization due to scattering over the im¬ 
pact polarization by the proton beam. The increase of tangential 
polarization in the FI model with the beam flux shows the com¬ 
pletely negligible effect of the impact polarization in comparison 
to other effects; the variation of the polarization is driven mainly 
by shifts of the optical depth scale due to higher ionization de¬ 
gree, which can lead to an increase of the resonance scattering 
effect on polarization. Another effect is the increase of radia¬ 
tion intensity due to non-thermal excitation in the lower chromo¬ 
sphere layers and an increase of the density of background col¬ 
liders. All these effects compete in the formation of these small 
fractional polarizations. The impact polarization contributes by 
only a small fraction to the total polarization degree in all the 
studied cases, as was expected fro m Fig. |71 

Contrary to the conclusions of IVogt et al.l (1200 ll) . the vari¬ 
ation of spectral index S does not affect the emergent polariza¬ 
tion degree much. This is because of the negligible impact polar¬ 
ization effect and because the number of low energy protons in 
the Ha formation layers is similar across the 6 values - the ini¬ 
tial difference in energy distribution of lower energy protons is 
quickly lost at depths ~ E^. An effect that could play a role is 
the variation of E^ energy. The energy distribution of the beam of 
protons is sensitive to the chromospheric depth parametrized by 
En (see Sect.|2|l. One expects that a fine tuning of E^ with the en¬ 
ergy necessary to reach the Ha-forming layers can lead to a more 
significant effect of impact polarization. We have calculated the 
dependence of the line center polarization on E^ for a wide range 
of cut-off energies. The resulfs can be found in Figs.fTOlandfTTl 
It is easy to find the effect of impact polarization that emerges at 
cut-off energies comparable to the energy necessary to reach the 
Ha-forming layer. This effect is seen as a decrease in the total 
polarization degree (i.e., increase of radial atomic polarization 
due to impacts that has an opposite orientation than the tangen¬ 
tial polarization formed preferentially by scattering). In the FI 
model, the thermal case of 0.01 % polarization is increased as 
protons have enough energy to change the structure of the atmo¬ 
sphere at Ha formation levels. The effect of impact polarization 
is seen at about 150 keV when the density of protons is high 
enough to decrease the dominant tangential polarization. Further 
increase of Ec again leads to a negligible number of protons with 
small energies at the Ha formation layer. Similar effects can be 
seen in the VAL F case, where the resonance polarization scat¬ 
tering effect is much stronger. The total number of thermal per- 
turbers is always lower than in the FI case, but high enough to 
decrease the thermal atmosphere polarization of 0.4 % even for 
a small energy cut-off. The impact polarization is best seen at 
energies around 200 keV; in this cooler model the total column 
mass above thq. = 1 is higher than in the FI case. In any case the 
emergent linear polarization is low, below the stated 5 % value, 
and it is even possible to expect a decrease of the total linear 
polarization due to proton impacts. 


9. Conclusions 

The purpose of this paper was to verify whether the anisotropic 
excitation of the n- 3 level of hydrogen caused by proton beams 
can lead to Ha line polarization and what degree of polarization 
can be expected. In our considerations, we have chosen a unidi- 
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Fig. 10. Ha line center polarization as a function of a lower 
energy cut-off Fc at the injection site. The atmosphere model 
FI, beam flux So = 10" ergcm^^ s“', S - 5. The E^ value has 
been varied in the interval 85-300 keV. 
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Fig. 11. Same as Fig. [TOl here for the chromosphere model 
VAL F and the E^ energy range 85-900 keV. 



rectional vertical beam with an initial power-law energy distri¬ 
bution at the chromosphere injection site, which is not deflected 
by collisions with constituents of the chromosphere although 
it is decelerated. For this extremely anisotropic beam we have 
calculated new chromosphere conditions with a wide range of 
beam parameters. For all the calculated models we have found 
the polarization degree well below the values reported in the 
past. Furthermore, the theoretical polarization degree is mainly 
caused by resonance scattering. 

The absen ce of impact polariz ation is consistent with mea¬ 
surements of iBianda et al.l (l2005h . although their conclusions 
about the beam isotropization are not necessary for explana¬ 
tion of the missing polarization. In fact, a sufficient number of 
unidirectional low-energy protons at the Ha-forming layer can¬ 
not be achieved for beams with initial power-law distributions. 
Although the proton beam may significantly affect the line in¬ 
tensity profile, impact atomic polarization is destroyed by colli- 
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sions with background electrons and protons and by the strong 
radiation held. The density of background perturbers is actu¬ 
ally higher than that calculated in previous works. In addition, 
the beam energy distribution at the line formation layers cannot 
be approximated by a power-law distribution and the effect of 
a steep increase of fractional polarization with spectral index 6 
is not observed. Therefore we believe that the polarization oc¬ 
casionally measured (i.e., by some authors) has another source 
than impacts by proton beams. 

Low impact polarization of the Ha line by proton beams in 
solar flares does not seem to be a good candidate for straightfor¬ 
ward proton beam diagnostics, at least if most of the usual initial 
distributions of proton beam energies are considered. It is possi¬ 
ble to imagine different initial energy distributions of the same 
energy flux with very narrow band of energies, which preferen¬ 
tially excite the hydrogen at low energies and in the Ha-forming 
layers. However, the question remains if such a distribution is 
physically possible and if it can lead to measurable impact po¬ 
larization. 

Recent sp ectropolarimetric measurements of an M6.3 flare 
performed bv IXu et ^ (l2005bl) indicate the presence of both ra¬ 
dial and tangential polarization of the Ha line and radial polar¬ 
ization of the Hy6 line at the edges of flare kernels. The process 
that is be lieved to play a significan t role is the effect of return 
currents (iKarlickv & Henou^r2002h . Investigation of this phe¬ 
nomenon remains the subject of further observational and theo¬ 
retical research. 
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